Functional renormalisation group approach to far-from-equilibrium quantum field dynamics 
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Dynamic equations for quantum fields far from equilibrium are derived by use of functional renormalisation 
group techniques. The obtained equations are non-perturbative and lead substantially beyond mean-field and 
quantum Boltzmann type approximations. The approach is based on a regularised version of the generating 
functional for correlation functions where times greater than a chosen cutoff time are suppressed. As a central 
result, a time evolution equation for the non-equilibrium effective action is derived, and the time evolution of the 
Green functions is computed within a vertex expansion. It is shown that this agrees with the dynamics derived 
from the 1 /A/'-expansion of the two-particle irreducible effective action. 
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Introduction. Far-from-equilibrium quantum field dynam- 
ics is one of the most challenging issues both in experimental 
and theoretical physics to date. Experiments exhibiting quan- 
tum statistical effects in the time evolution of many-body sys- 
tems are extremely demanding. In particular, the preparation 
of ultracold atomic Bose and Fermi gases in various trapping 
environments allows to precisely study quantum many-body 
dynamics of strongly correlated systems, see, e.g., Refs. IjI]. 
In recent years, the field has attracted researchers from a va- 
riety of disciplines, ranging from condensed-matter to high- 
energy particle physics and cosmology. Nonequilibrium field 
theory to date is dominated by semi-classical mean-field ap- 
proaches which are in general only valid for weak interactions 
or large occupation numbers. For strongly correlated quantum 
systems methods are available predominantly for systems in 
one spatial dimension and include the Density Matrix Renor- 
malisation Group methods, see e.g. 12], as well as techniques 
for exactly solvable models |3]. Non-perturbative approxima- 
tions of the two-particle irreducible (2PI) effective action |3l 
have been intensively studied and applied to nonequilibrium 
dynamics flHSHSIlS U 1 and are applicable also in more 
than one dimension. In field theory, they also provide a way 
to study strongly correlated fermions beyond mean-field and 
Boltzmann approximations |8]. Like these, the results pre- 
sented here are expected to be of high relevance, e.g., for the 
description of ultracold degenerate Fermi gases close to the 
BEC-BCS crossover 1 12]. 

In this paper we derive dynamic equations for quantum 
fields far from equilibrium. As a central result, we derive a 
new time evolution equation for the non-equilibrium effec- 
tive action by use of functional renormalisation groiip (RG) 
techniques, cf. Refs. (HI H HI [l6t], as well as lO for 
non-equilibrium applications. This exact and closed evolu- 
tion equation allows for non-perturbative approximations that 
lead substantially beyond mean-field and quantum Boltzmann 
type approaches: it can be rewritten, in a closed form, as a hi- 
erarchy of dynamical equations for Green functions, and this 
hierarchy admits truncations that neither explicitly nor implic- 
itly rely on small bare couplings or close-to-equilibrium evo- 
lutions. Moreover, the hierarchy of equations is technically 
very close to 2PI and Dyson-Schwinger type evolutions, see 
lll6ll . as well as to evolution equations for the effective action 



as derived in |18|] . This has the great benefit that results from 
either method can be used as an input as well as for relia- 
bility checks of the respective truncation schemes. In turn, in 
particular the different resummation schemes invoked in these 
approaches, as well as the differing dependences on Green 
functions, allow for a systematic analysis of the mechanisms 
of non-equilibrium physics. Here we apply the above setting 
within a non-perturbative vertex expansion scheme as well as 
an s-channel approximation. This truncation turns out to cor- 
respond to an expansion in inverse powers of the number of 
field components M. The dynamic equations derived agree to 
next- to-leading order with those obtained from a expan- 
sion of the 2PI effective action. 

Functional renormalisation group approach. For a given 
initial-state density matrix pnito), the renormalised finite 
quantum generating functional for time-dependent n-point 
correlation functions, 



(1) 



(summation over double indices is implied) carries all the 
information of the quantum many-body evolution at times 
greater than the initial time to ■ The Heisenberg field operators 
^a{x), a = 1, ...,Af, are assumed to obey equal-time bosonic 
commutation relations. In Eq. ([T]|, Tc denotes time-ordering 
along the Schwinger-Keldysh closed time path (CTP) C lead- 
ing from to along the real-time axis to infinity and back to to, 
with J^Q = dxQ J d'^x. All connected Greens functions 
will be time-ordered along C. 
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FIG. 1: (color online) The closed time path C(r) terminating at the 
time given by the parameter r. At later times, field fluctuations 
summed over in the generating functional do not contribute to Green 
functions the maximum time argument of which is r. 

The key idea of our approach to dynamics is to first consider 
the generating functional for Green functions where all times 
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are smaller than a maximum time r. This implies a time path 
C(r) which is closed at t = r, see Fig.[Tl and we are led to the 
generating functional = ^c(t) with the source term 



(2) 



At r = to, this results in a trivial Ztg where all information 
is stored in the initial density matrix poito)- From this initial 
condition Zr can be computed by means of the time evolution 
drZr for all times t > to. 

We will derive this evolution by using functional RG ideas. 
To that end we note that Z^ can be defined in terms of the 
full generating functional Z^o in ([T]) by suppressing the prop- 
agation for times greater than r. This suppression is achieved 
by 



functional (O is equivalent to a generating functional with a 
closed time path C(t) leading from to to r and back to to- 
Note that the CTP automatically arranges for the normalisa- 
tion of Zt-. We conclude that, by construction, the sharp cutoff 
(|4]l is a physical one in the sense that it corresponds to inte- 
grating out all fluctuations being relevant for the evolution up 
to a particular time. 

The restriction of the CTP to times to < t < t implies 
that the differential equation for Zr describing the flow of 
the generating functional, and therefore that of the correla- 
tion functions, encodes the full time evolution of the sys- 
tem. Analogously, the time evolution of connected correla- 
tion functions is derived from that of the Schwinger functional 
Wt = —ilnZr- It is more convenient, however, to work with 
the effective action. 



(5) 



where the function is chosen such that it suppresses the 
fields, i.e., 6/6 Ja, for all times t > t. This requirement does 
not fix Rr in a unique way, and a simple choice is provided 
by 



iRr,abix,y) 







for xo = yo > T,x 
otherwise 



y, a 



(4) 

see Fig. 121 We note that for t oo, we recover the full gener- 
ating functional. Zoo = Z, while for t — to, any evolution to 
times t > to is suppressed. In the latter case all rt-point func- 
tions derived from Zr are defined at to only and reduce to the 
free classical ones, or any other physical boundary condition. 
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FIG. 2: (Color online) The cutoff function Rr,ab(x,y) in the time 
plane {xo,?/o} = t,t' > to- The function vanishes every- 

where except for t = t' > r where it tends to infinity and therefore 
implies a suppression of all fluctuations in the generating functional 
at times greater than r. 

We emphasise that the cutoff Rr in (O suppresses any time 
evolution at times greater than t. Correlation functions de- 
rived from Zr vanish as soon as at least one of their time ar- 
guments is larger than r. Hence, the regularised generating 



Here, space-time arguments are suppressed, and (pa{x) — 
6Wr/6Ja{x)\j=o is the classical field expectation value. 
From Eqs. (O and (|5]l we derive the Functional RG or flow 
equation for the r-dependent effective action. 



drTr ^ - 



df Rr,ab 



(6) 



where F^"^ = 5'^Tr/{6(l))'^. Again, space-time arguments are 
suppressed which appear in analogy to the field indices a, b, 
see e.g. Ref. Hi]. Eq. a represents our central result. It is 
analogous to functional flow equations used extensively with 
regulators in momentum and/or frequency space to describe 
strongly correlated systems near equilibrium 113, 14, f5, 3. 
Its homogenous part relates to standard r-dependent renor- 
malisation [16J, and has been studied eg. in Sim. We 
close the derivation of the evolution equation with some re- 
marks: Eq. (|6]l is nothing but an infinitesimal closed time loop 
at r, which is by itself finite and requires no further renormali- 
sation. However, in particular in higher dimensions additional 
regulators in the spatial or momentum domains may be advan- 
tageous for facilitating renormalisation and thus the practical 
application of the approach, see |jT3, 14, T^. We also em- 
phasise that the evolution equation (|6]l, even though close in 
spirit and construction to standard functional RG equations, is 
conceptually different. It entails a physical time evolution as 
opposed to integrating out degrees of freedom. Nevertheless, 
more general regulators Rr may be advantageous in other 
cases lll6ll . e.g., when a cutoff is set in the relative time (t—t'-) 
direction, see Fig.|2] Such a case would correspond to setting 
a cutoff along the frequency axis [17] which is conceptually 
closer to the standard RG approach with cutoffs in the mo- 
mentum domain. A more detailed discussion of our approach 
is deferred to 1 19]. 

Flow equations for correlation functions. To obtain a prac- 
tically solvable set of dynamic equations, we derive the flow 

in) 

equation for the proper n-point Green function Tr by taking 
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the nth field derivative of Eq. Fig.|3]shows a diagrammatic 
representation of the resulting equations for the r-dependent 
proper two- and four-point functions. 



(1) _ ' 



a.r,a.[*=0] = -i 




b' V 



FIG. 3: (color online) Diagrammatic representation of the general 
flow equations for r$-^^ [(^], T^^Vl- and rl^l^^^[(fi = 0], for a 0*- 
theory. Open circles with a r denote drRr.ab, solid lines with (blue) 
filled circles are r- and, in general, (/)-dependent two-point functions 
GT,ab = jfri^' + Rt]~^ ■ All Other filled circles denote proper field- 
dependent n-vertices T^"^^^^, n = 3, 4, 6. P implies a sum corre- 
sponding to all permutations of its arguments. 



We would like to emphasise that proper 7i-point Green func- 
tions in general do not vanish for times greater than t as op- 
posed to the correlation functions derived from Zr- However, 
in contrast to usual renormalisation-group flows in the mo- 
mentum domain, causality prevents the influence of proper 
Green functions for times greater than r on the dynamics up 
to time r. 

To be more specific, we consider in the following the spe- 
cial case of an TV^-component scalar 0"* theory defined by the 
classical action 



(7) 



where (^^ = ^a'fa, and where the inverse free classi- 
cal propagator involves one or more time derivatives, e.g., 
'^0^6(2^' y) = '5c(a; " y)\~'^'ib^xo + i^^iB(a;)(5afc] for a non- 
relativistic trapped Bose gas, with being the Pauli 2-matrix 
and IIyb,{x) — —/\\/2m + V{x) the one-body Hamiltonian. 

Our goal is to derive the full time-evolution of F*^") = 
r^-*, in particular, of the connected two-point function G = 
j[p(2)]-i = i[rM'']~^ which contains all information about 
the normal and anomalous one-body density matrices (cf., 
e.g. Ref. iid). As a consequence of the above mentioned ef- 
fective cut off of the CTP at times greater than t, it will be 
sufficient, for the time evolution up to i = t, to determine the 
functions Tr and thus the propagator 



We restrict ourselves to the case 

.(3) 



0, such that the ac- 

(1) 



tion (|7]) implies that Fi- = 0, and thus the flow of Fr van- 
ishes. Moreover, the equation for the proper two-point func- 
tion involves, on the right hand side, only the term containing 



T,ab ^ 2 



T,abcd 



{Gr [dTRT]Gr)dc , 



(9) 



see also Fig. |3] The term in parentheses stands for the regu- 
larised line. We supplement Eq. ^ with the flow equation for 
T'"r\ which, for (/)a = 0, is depicted in Fig. [3] This system of 
equations is still exact. For practical computations it needs to 
be closed which can be achieved by truncation or by supple- 
menting it with equations for one or more higher n-vertices 
truncated at some higher order Here we truncate by neglect- 
ing, in the equation for Fr^-*, the term involving F^^^ 



^■r^ T,abcd — g J \^ T,abef'~' -r,} 9^ r,cdgh j 



X {Gr[drRr]Gr)he + P(a,5,C,d). (10) 

P implies a sum corresponding to all permutations of its ar- 
guments. In this way we obtain a closed set of integro- 
differential equations for the proper functions up to fourth or- 
der As we will show in the following, they allow to derive, for 
a particular cutoff time r, a set of dynamic equations describ- 
ing the time evolution of the two- and four-point functions up 
to time t — T. We emphasise that the only approximation here 
is the neglection of the six-point vertex, see Fig. [3] 

Dynamic equations. For the sharp temporal cutoff Rr cho- 
sen here the flow equations can be analytically integrated over 
T. As pointed out above, our cutoff implies the connected two- 
point function to vanish at times greater than r, i.e., it can be 
written as 



F(2 



ab 



ta)e{T-tb), 



(11) 



(8) 



where 6{t) evaluates to for r < and to 1 elsewhere, and 
where ta is the time argument corresponding to the field in- 
dex a, etc. Hence, the precise way in which the cutoff Rr^ab 
diverges Sitta = h > t is chosen such that Fr +Rt is the in- 
verse of ~iGr for all times ta, tb, see Eq. (O. Using Eq. (fTTT i 
one finds that 

{Gr[drRr]Gr)ab = -lGr.abdr[eiT - ta) 0(t - tb)]. (12) 

Note that we have not used the specific choice (|4|l for deriving 
(fT2b but simply the property ( fTTT ). that is the suppression of 
any propagation for times t > t. This independence of the 
specific choice of Rr leading to (fTTI) is the generic feature of 
the present approach. After inserting Eq. (fT2] | into Eqs. ^ 
and ( [Tol l, we can integrate over r and obtain, after some alge- 
bra, the integral equations determining the flow of the proper 
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functions from to to some final time t, 

t \ ft 
to 2 



.(2) 
T.ah 



.(4) 
r.abcd 



tafi 
t 

tafi 



(13) 



r 



(4) 



X r 



rl,,,cdahGr^uM + (a ^ c) + (a ^ d). (14) 



Double indices imply sums over field components, spatial in- 
tegrals and time integrations over the CTP C, from to t and 
back to io. We furthermore introduced 



Tab = max{ta,tf,}, 
Tabcd = max{ia, tf,, ic, id}- 



(15) 



The brackets denote terms with the respective indices 
swapped. 

From Eqs. (fTsT l and (fT4l i it is clear that for the two- and 
four-point functions to be defined at r = t we need to specify 
initial functions air — t^. We point out that, within the trun- 
cation scheme chosen above, we can insert any set of proper 
two- and four-point functions defined in their time arguments 
at and only at to. as long as we set all n-vertices for n = 1, 3, 
and n > 4 to vanish. Our scheme corresponds to a Gaus- 
sian initial density matrix pnito) since the four-point func- 
tion has an influence on poit) only for t > to. Here, we 
choose the respective classical proper functions defined by S 
in Eq. (|7]i. Hence, the initial two- and four-point functions 
entering Eqs. ( fT3] l and (fT4l l read 



r 



(2) 
tQ,ab 



(2) 



O.ab' 



(16) 



^Itbcd = ^ilL = -i^g/mSabScd + SacSbd + SadSbc) 

X dc{xa - Xb)Sc{xb - Xc)5c{xc - Xd). (17) 

In order to arrive at a set of dynamic differential equations, we 
finally rewrite Eq. (fTsl l as 

^'Ga.acGT,b,cb = «<5c,afc " 7^ j ^Ta^.adce Gra^.ed Gr^^,cb, 
^ Jto,C 

(18) 

with 5c,ab = 5ab5c{Xa " Xb) = SabSc{ta - 4)<5(Xa - Xf,). 

Eq. ( fTSl) is the dynamic (Dyson-Schwinger) equation for the 

(2) — 1 

connected two-point function Gr^^^ab- Since T).^' = iG^ 
represents a differential operator, the solution of Eq. ( fTSl ) fi- 
nally requires another set of boundary conditions to be spec- 
ified, depending on the form of the differential operator. In 
our case, this is the initial two-point function, fixed by the 
one-body density matrices as well as Bose statistics 19(]. We 
finally point out that we have set t — Tcb since the flow of 
Gr.cb Stops at the maximum of the time arguments tb, t^- This 
can be proven from the structure of Eqs. ( fT4b . ( fTSl ) but is more 




(tf, tg = (■„... mm(r',0 ) 




(t^ = t„... max(r,„g ) 

FIG. 4: (color online) The upper equation is the s-channel projec- 
tion of Eq. l ll4b . The second equation defines the resummed ver- 
tex appearing in the Dyson-Schwinger equation derived from the 
NLO l/Af approximation of the 2PI effective action | 5]. The two 
definitions are identical in every order of a perturbative expansion 
(see text). Dashed lines denote the s-channel part of the bare ver- 
tex rfj,^;,, see text. All other symbols correspond to those in Fig. [3] 
Letters on internal lines indicate summation over field indices and in- 
tegration over space and time (along the CTP from to to t and back). 
The integration intervals are given in parentheses. 



easily seen from the definition (O: Once the hard cutoff t has 
passed the largest time appearing in a (connected) correlation 
function, the flow stops since the forward and backward parts 
of the CTP over greater times cancel identically in the func- 
tional integral. 

Let us assume that t = ta denotes the present time, at which 
Eq. (fTSl ) determines the further propagation of Gt^,ab for tb < 
ta (for tb > ta the solution is then fixed by symmetry). We 
point out that all time arguments of the functions occurring on 
the right-hand sides of Eqs. ( fT3] l, ( fT4b . and therefore all cutoff 
times Tab and T^fgh are evaluated at times t' < t. Hence, 
in accordance with causality, Eqs. ( fT4l l and ( fTSl l for a given 
initial Gto.abitoito), can be solved iteratively. 

This concludes the derivation of our main results, a closed 
set of dynamic equations for the two-point correlation func- 
tion as obtained from a functional RG approach with a cutoff 
in real time. In the remainder of this article we will concen- 
trate on rederiving the intensively used dynamic equations ob- 
tained from the 2PI effective action in NLO of an expansion 
in inverse powers of the number of field components J\f and 
comment on the potential of our approach beyond this approx- 
imation. 

From RG to 2PI next-to-leading order (NLO) 1/Af. The 
mean-field approximation, i.e., the well-known Hartree-Fock 
equations are obtained by neglecting the flow of the four-point 



according to Eq. ( fT4] i. F 



(4) 



As a consequence, the 



flow parameter of Gr,ab is, for the chosen cutoff, fixed to r = 
Tab and can thus be neglected. Cf. Ref. |0] for the mean-field 
equations in the notation used here. 

As a first step beyond mean-field we consider the trunca- 
tion in which the s-channel scattering diagram is included be- 
yond the mean-field limit in the loop integral on the right- 
hand side of Eq. ( fT4] i. The result of this section will be 
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that the obtained equations are equivalent to the 2PI equa- 
tions in NLO of a 1/M expansion [iSJ. This truncation cor- 
responds to keeping only one channel of each, the classical 
vertex and the one-loop integral term. The four-vertex then 
only depends on two space-time variables and field indices, 



^ t,acbd ^i,ab^o-c^bd, and enters, in Eq. ([T8]l, as a cutoff- 
dependent self energy Sr^t^ab = ^^r^^.ab'^T^t.ab- We can al- 
ways write t — min(t,i'), with t = t' = Tab, for the cutoff 
parameter in r|^^^, see Fig. |4l upper line. Since the param- 
eter in the (green) vertices on the right hand side is also the 
minimum of the maxima of integration times in the adjacent 
loop and the respective external time t or t' , one can iterate 
the integral equation in order to obtain a perturbative series 
of bubble-chain diagrams consisting only of classical vertices 
and full, cutoff-dependent propagators. This procedure pro- 
vides us with a proof that E^^j, is identical to the NLO 
2PI 1/M self energy obtained from the 2PI effective action, 
cf. Ref. 15s] and Fig.|4] lower line. 

This proof which exploits the perturbative expansion is 
lengthy and will be given elsewhere. Here we show that the 
above identity can be inferred in a comparatively easy way 
from the topology of the different terms in the flow equations 
for the two-, four- and six-point functions: Consider the un- 
truncated set of equations as displayed in Fig. |3] First, non- 
s-channel contributions do not generate bubble-chains of the 
form shown in Fig.|4] Second, is one-particle irreducible. 
Its contribution to the flow of F^^^ does not give rise to bubble- 
chains, even if inserted recursively into the first diagram on 
the right-hand side of the flow equation for F^^^ In turn, 
by dropping the second diagram with the six-point function 
and using the s-channel truncation, the iterated flow equation 
generates only bubble-chain diagrams with full propagators as 
lines. Hence, a r-integration of this set of flow equations leads 
to dynamic equations which include all bubble-chain contri- 
butions and are therefore equal to those obtained to NLO in a 
1/N expansion of the 2PI effective action |5]. We emphasise 
that the above topological arguments are generally valid when 
comparing resummation schemes inherent in RG equations of 
the type of Eq. Q, with those obtained from 2PI effective ac- 
tions. This applies , e.g ., to equilibrium flows iHlllIii] and 
thermal flows il5ll22ll . For a comparison with 2PI results see 



Ref. 1 22 1, for the interrelation of 2PI methods and RG flows 
Ref. lll6ll . In this context we would also like to remark that 
our approach, by constructions, obeys Goldstone's theorem as 
it is solely based on the IPI effective action and its derivatives, 
the proper n-point functions. Hence, spontaneous symmetry 
breaking always goes together with gapless excitations in the 
proper 2-point functions. For a resolution of this matter in the 
2PI approach see e.g. Refs. ll6. 231. 

Equilibration of a strongly interacting Bose gas. We close 
our paper with an illustration of the non-perturbative nature 
of the dynamic equations obtained in the s-channel approx- 
imation, which are, as shown, equivalent to the 2PI equa- 
tions in NLO llN approximation. For this we face the long- 
time evolution of a weakly interacting non-relativistic Bose 
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FIG. 5: Time evolution of momentum-mode occupation numbers 
n{t; p) /niL of a uniform ID Bose gas, normalised by the total num- 
ber niL of atoms in the box. The initial far-from-equilibrium state 
corresponds to a Gaussian momentum distribution. Panel (a) de- 
scribes the evolution of a weakly interacting gas, 7 — 1.5 • 10~^, 
panel (b) that of a strongly interacting gas with 7 = 15. 



gas obtained in Ref. flTl, for which 1/7V = 1/2, with that 
of a strongly interacting one. These results illustrate, what 
is supported by the non-perturbative approach presented here 
as well as by benchmark tests for special-case systems ifj^ll . 
viz. that the s-channel or 2PI 1/JV approximation is applica- 
ble for strong coupling also for small J\f where an expansion 
in powers of 1 / J\f would seem questionable. The uniform gas 
in one spatial dimension starts from a non-equilibrium state 
which is assumed to be fully described by a non-vanishing 
two-point function Gaa(ioj*o;p) = n{tQ;p) + 1/2, a = 1,2. 
The initial momentum distribution of atoms in the gas is 
n{to;p) = (ni0rcr) exp(-p^/cr^), with o- = 1.3- 10^ m~^. 
The computations were done as described in Ref. [11], on a 
grid with Ng = 64 points spaced by Ug = 1.33/im. Fig.|5] 
shows the evolution for a (a) weakly interacting sodium gas, 
with line density ?ii = 10^ atoms/m and interaction strength 
g — fi^jni/m, 7 = 1.5 x 10^'^ and (b) a strongly interacting 
one with 7 = 15, rti = 10^ /m. The final momentum distri- 
butions closely resemble Bose-Einstein, i.e., thermal distribu- 
tions n{t;p) oc l/[cxp{f3[u!p — /i]) — l], with the dispersion 
Ljp extracted from the numerical data 01 III . 

Conclusions. Using functional RG techniques we have de- 
rived a theory of far-from-equilibrium quantum field dynam- 
ics. The time evolution of the system is generated through the 
flow with a cutoff time introduced in the generating functional 
for correlation functions. In a truncation of the flow equations 
we recover the dynamic equations as known from the 1 /J\f- 
expansion of the 2PI effective action. These results lead us 
to the following conclusions: Although sub-leading terms in 
the expansion of the 2PI effective action are suppressed with 
additional powers of 1 /J\f they contain bare couplings, and, 
despite their resummation, the approximation seems formally 
questionable for strong couplings if g/Af is large. Presently, 
an intense discussion focuses on the question to what extent 
the NLO 1/JV approximation is applicable for large interac- 
tion strengths, drawing from benchmark tests for special-case 
systems [24]. Our results provide analytical means to find 
the conditions for the validity of this approximation since no 
perturbative ordering in bare couplings is involved, and the 
remaining truncation in the Green functions can be tested for 
self-consistency. 
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